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ABSTRACT 

Extragalactic sources from the IRAS Faint Source Catalog (FSC) which have 
the optically faintest magnitudes {E > 18) were selected by spatial coincidence 
with a source in the FIRST radio survey, and 28 of these sources have been 
observed with the Infrared Spectrograph on Spitzer (IRS). While an infrared 
source is always detected with the IRS at the FIRST position, only ~ 50% of 
the infrared sources are real FSC detections, as estimated from the number of 
sources for which the fu{25 fim) determined with the IRS is fainter than the 
sensitivity limit for the FSC. Sources have 0.12 < z < 1.0 and luminosities (ergs 
s~^) 43.3 < log[z/Li,(5.5 yum)] < 46.7, encompassing the range from local ULIRGs 
to the most luminous sources discovered by Spitzer at z ~ 2. Detectable PAH 
features are found in 15 of the sources (54%), and measurable silicate absorption 
is found in 19 sources (68%); both PAH emission and silicate absorption are 
present in 11 sources. PAH luminosities are used to determine the starburst 
fraction of bolometric luminosity, and model predictions for a dusty torus are used 
to determine the AGN fraction of luminosity in all sources based on z^L,^(5.5 /im). 
Approximately half of the sources have luminosity dominated by an AGN and 
approximately half by a starburst. The ratio of infrared to radio flux, defined as q 
= log[/,y(25 /im)//i,(1.4 GHz)], does not distinguish between AGN and starburst 
for these sources. 

Subject headings: infrared: galaxies — galaxies: starburst — galaxies: AGN 

1. Introduction 

As part of efforts to determine the nature of faint infrared sources, hundreds of spectra 
have been obtained with the Infrared Spectrograph on Spitzer (IRS; Houck et al. 2004) 
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of optically faint sources discovered at 24 fim with surveys using the Multiband Imaging 
Photometer (MIPS; Rieke et al. 2004). The most important result is that these faint sources 
having f,^(24 nm) ~ 1 m Jy are generally at high redshift, typ ically z ~ 2 (IHouck et al.ll2005 



Weedman eraPboOGal lbl: IVan et alj[2007l : ISrand et allboO?! ). Most show the 9.7 ;um silicate 
absorption feature, but some have PAH emission features and some have only a continuum 
with no detectable features. These results imply that the sources are optically faint because 
they are both distant and obscured by dust. 

The absorb ed sources and the featureless, power-la w sources have been interpreted as 
obscured AGN ( ISajina et al.l 120071 : iPoUetta et al.l 120071 ). This interpretation implies that 
the surface density of optically obscured AGN at high redshifts exceeds that for classical, 
opticall y discoverable AGN by a factor of two to three in infrared su rveys to /i.(24 ^m) < 



1 mJy jWeedman et al.ll2006cl : IPoUetta et aDl2007l: iFiore et al.1 120071'). Extre melv luminous 



starbursts have also been found at z ~ 2 JParrah et al.l 



2008 



Pope et al.ll2008l ) which appear 



similar to local starbursts in all characteristics except luminosity. 

To understand these high redshift, dusty sources and to interpret the Spitzer-discoYeied 
sources within scenarios for evolution in the universe, it is essential to understand how 
the sources of luminosity are distributed between AGN and starbursts among the overall 
population of sources. This requires comparison with closer, brighter examples for which 
AGN and starburst diagnostics within the IRS spectra can be compared to information from 
other wavelengths. So far, the compar i son samples which show the most sira i lar spectral 
shap es are the ULIRGS (jllao et al. 2007; Imanishi et al. 2007; Farrah et al. 2007; Desai et al. 



20071 ) , but the ULIRGs have lower luminosities and smaller redshifts than the sources at z ~ 
2. The ULIRGS observed with the IRS were chosen primarily on the basis of 60 /im fluxes 
from the Infrared Astronomical Satellite (IRAS) and were not selected based on optical 
characteristics. 

In this paper, we describe a sample of optically faint, infrared-bright galaxies which spans 
luminosities and redshifts between the ULIRG and high redshift Spitzer IRS samples. This 
new sample ari ses from selecting th e optically faintest sources listed in the IRAS Faint Source 
Catalog (FSC) (iMoshir et al.lll990l ). These sources have E > 18 mag (the photographic E 
magnitude band is cent ered at 0.66 fim and approximately corresponds to R magnitude 
( jHumphreys et al.lll99ll )). 



(Becker et al. 


1997; 


White et al. 


1997) 



detection was used to increase confidence that an FSC source is real (as opposed to noise or 
cirrus) and to provide an accurate position for selection of the optical counterpart. 
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New IRS spectra of 24 sources and 4 spectra of archival sources are presented. This 
sample of 28 sources is compared to previous samples of bright IRAS sources for which 
Spitzer IRS spectroscopy is available, and the validity of the FSC identifications is discussed. 
Spectral parameters are used to estimate starburst and AGN contributions to the luminosity. 



2. Sample Selection and Data Analysis 

2.1. Selection of New Sample of Infrared Bright, Optically Faint Sources 

To select a sample of optically faint infrared sources from the IRAS FSC, we first 
cross-correlated the IRAS FSC sources with the FIRST catalog, using FSC sources listed 
as real detections rather than upper limits at both 25 ^m and 60 ^m (sources with quality 
flags 2 or 3). 2,310 sources were associated to within 3 times the FSC one sigma positional 
uncertainty (ranging from 15 "to 90 "for individual sources). Images from the Palomar Second 
Digitized Sky Survey (DSS2) were then examined for the 2,310 fields to determine an optical 
identification at the radio source position. 

Of the 2,310 sources, 1,944 can be identified with bright or medium brightness galaxies 
(photographic E < 16 mag), and 225 are fainter galaxies (16 mag > E < 21 mag). The 
remaining 141 sources are identified with bright Galactic stars or have no optically identifiable 
galaxy; the latter could be either spurious sources or cirrus sources. 

The best available optical photometry for the entire FSC sample is the Minnesota APS 



catalog (jCabanela et al.ll2003l ). which gives photographic O and E magnitudes. The final 
sample we select is defined by having E > 18. While actual magnitudes are uncertain to 
a few-tenths, we can have confidence that the APS magnitudes are homogeneous and allow 
location of the optically faintest sources. The selection for E > 18 yields 31 sources, which 
is about 1% of all sources common between FSC and FIRST. 



2.2. IRS Observations and Analysis 

We have obtained IRS low-resolution spectra for 24 of these 31 sources, and an additional 
4 have been observed in other Spitzer programs. The sample of objects and the details of 
the Spitzer IR^ observations are in Table 1. Our new observations were made with the 
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IRS Short Low module in orders 1 and 2 (S Ll and SL2) a nd wi th the Long Low module in 
orders 1 and 2 (LLl and LL2), described in iHouck et al.l (120041 ) . These give low resolution 
spectral coverage from ~5 /im to ~35 /im . 

Sources were placed on the slit by using the IRS peakup mode with the red camera. 
Background subtraction for spectra was done using coadded background images that added 
both nod positions having the source in the other slit (i.e., both nods on the LLl slit when 
the source is in the LL2 slit produce LLl spectra of background only). This subtraction 
is possible when both SL orders or both LL orders are observed with the same integration 
times. Independent extractions of spectra at the two nod positions were compared to reject 
any highly outlying pixels in either spectrum, and a final mean spectrum was produced. 

Starting with vl5 of the S SC basic calibrated data (BCD), spectra were extracted using 
the SMART analysis package (iHigdon et al.ll2004l ). Extractions were done with an average 
width of 4 pixels (perpendicular to dispersion; width varies with wavelength because of 
varying spatial resolution). Calibrations applied to BCD products use extractions of 8 pixel 
width for the calibrating stars, so we correct fluxes by using Markarian 231 as a standard 
source and extracting it with both 8 pixel and 4 pixel windows. Final spectra were boxcar- 
smoothed to the approximate resolution of the different IRS modules (0.2 /im for SLl and 
SL2, 0.3 yum for LL2, and 0.4 /im for LLl). All spectra are illustrated in Figures 1-4, ordered 
by luminosity. 

Measured continuum parameters and feature strengths for the spectra are in Table 2. 
Fluxes and equivalent widths (EWs) for the 6.2 /im and 11.3 /xm PAH features are measured 
using the single Gaussian line fit routine within SMART. For the 6.2 fim PAH line, the fitting 
is between 5.5 /im and 6.9 /im ; for the 11.3 /im line, fitting is between 10.4 /im and 12.2 /im . 



As discussed by iBrandl et al.l (120061 ) and IPope et all (120081 ). fluxes and EWs of PAH 
features depend on the underlying continuum level which is selected. This is not a straight- 
forward choice and is more uncertain for the blended 7.7/im feature than for the individual 
6.2 /im and 11.3 /im features. To minimize uncertainties based on assumptions for fitting the 
features, we do not list measurements for the 7.7/im line. 

The fitting procedure we use applies a linear fit for the continuum and a Gaussian 
fit for the line between the selected wavelengths. This l inear continuum fit is somewhat 
different from the spline fit used by iBrandl et al.l (120061 ) and yields a higher continuum 
level. To estimate uncertainties arising from the fitting procedure, we remeasured the Brandl 
starbursts with the procedure we use. Using the 17 starbursts having single components (so 
there is no uncertainty regarding w hat flux to inc l ude in the spectral extraction), the average 
EW of the 6.2 /im feature given by IBrandl et al.l (120061 ) is greater by 1.11 then the EWs we 
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measure. 



In Table 2, the parameter r^j describing the depth of the 9.7 /ir n sihcate absorptio n 
feature is measured using the continuum fitting technique described by lSpoon et ahl (120071 ). 
The average spectrum of the 19 sources showing sihcate absorption is compared to the 
spectrum of Markarian 231 in Figure 5. 

We note that Source 9 is a blazar (^Perry et ah I ll978h with observed IRS flux at 25 /im smaller 
by a factor of 30 than when observed by IRAS. The source met the criteria for definition of 
the sample and so was observed, but this source is not included in the luminosity plots. As a 
result of the unexpectedly low flux, the IRS spectrum is very noisy and shows no measurable 
features (Figure 4). PAH limits are not meaningful for this source compared to those for 
other sources. 



3. Discussion 

3.1. Reality of IRAS FSC Sources 

The optically faint sources chosen are often near the lower limit of /i/(25 fim) in the 
IRAS FSC; the mean FSC /i.(25 /xm) of the sources chosen is 150 mJy, but the typical one 
sigma noise for an FSC source is 65 mJy. At this low signal to noise, some of these sources 
may be spurious. They could appear in the FSC catalog either because of enhancing noise 
or because of confusing Galactic cirrus. 

To test the reality of the FSC sources, we determined what should be the IRAS /^(25 
/im) by applying a synthesized IRAS filter to the Spitzer IRS spectrum, using the synthesis 
tool in SMART. The resulting IRS /iy(25 /im) are listed in Table 2. In Figure 6, the IRAS 
/^(25 /im) and the synthetic IRS /,y(25 fim) are compared. 

Fig. 6 shows that the /;y(25 /im) from IRAS are always greater than the f^{25 fim) 
measured with the Spitzer IRS, with a median factor of 1.6 and extending to a factor of 6. 
The discrepancy increases with decreasing flux, so that the faintest sources show the largest 
ratio of IRAS to Spitzer flux. Various previous observations with comparable IRS signal 
to n oise indicate inter nal consistency of Spitzer MIPS and IRS fluxes at 24 /im to within 



5% (IHouck et al.l 120071 ). so IRS flux calibration uncertainty is much less than the differences 



observed with IRAS fluxes. 

One possible source of difference arises when the IRAS source is larger than the IRS 
slit so that so that s ome flux is omitted in the IRS spectrum. This is the case for many 



nearby starbursts in iBrandl et al.l (120061 ). Images of the FSC sources were obtained at 
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25 yum because acquisition of sources used the IRS red peakup camera, and all images show 
an unresolved source as evidenced by the presence of a diffraction ring. This indicates that 
extended sources cannot explain the differences between IRAS and IRS fluxes. 

We conclude that the differences between IRAS and IRS fluxes arise because many of 
these sources are near the detection limit for the IRAS FSC Fig. 6 shows that 14 of the 
IRAS FSC sources are actually fainter at 25 fim, as measured by IRS /jy(25 fim), than the 
typical one sigma sensitivity limit of 65 mJy quoted for the FSC. This implies that these 
faint FSC "detections" are spurious, perhaps accidental noise spikes. 

Therefore, we use the ratio IRAS /i/(25 /im)/IRS fu{25 jim) to estimate the reality of 
the IRAS FSC source. If this ratio is < 1.5, we assume the source is really detected in 
the IRAS FSC, because flux discrepancies of this amount are reasonable for the faint FSC 
sources. Also, this is the approximate flux ratio at which the IRAS FSC fluxes fall below the 
65 mJy sensitivity limit in Figure 6. By this criterion, 14 (50%) of the sources are real FSC 
detections (counting the blazar for which the flux discrepancy is attributed to variability). 

If half of the FSC infrared sources are not real, why did all of the FSC sources appear 
as Spitzer sources? We believe that the answer is because the requirement of a FIRST radio 
source at the IRAS FSC position was also a criterion for source selection. As shown below, 
Spitzer sensitivities are such that a Spitzer source should be detected for any radio source 
brighter than the FIRST limit, regardless of whether the source was actually detected in the 
IRAS FSC. 

The ratio of infrared to radio flux is described by the parameter q = log[/,y(25 /xm) / /iy(1.4 
GHz)]. The median q = 0.8 for faint Spitzer sources detected at both 24 /im and 1.4 GHz 



(lAppleton et al.l 120041 ). and this is a lso the q value expe cted from previously known radio 



infrared correlations for starbursts (jCondon et al.l 119821 ). The FIRST catalo gue includes 



sources to a point-source limit of 1 mJy for typical one sigma rms of 0.13 mJy (jWhite et al. 



19971 ). This limit means that a Spitzer source would have been found with /i/(25 jim) > 6 
mJy for any FIRST source with q > 0.8. Such a source would be detectable by the peakup 
pointing of Spitzer, and an IRS spectrum would be obtained, even if the FSC detection were 
not real. 

This result means that any FIRST source with an infrared to radio ratio exceeding 
the median expected value would be detected in our Spitzer pointings. The existence of a 
FIRST source makes it probable, therefore, that a detectable Spitzer source would be found 
even if the FSC source used for the initial identiflcation were not real. This explains why 
Spitzer sources were always found at the FIRST position even in those cases when the fiy{25 
fim) are too faint for the FSC detection to be real. 
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3.2. Characteristics of Spectra and Near-Infrared Luminosities 



Spectra of the FSC sample are shown in Figures 1 through 4, ordered by z/L,y(5.5 ^m), 
with lowest luminosity sources shown first. This ordering allows some overall spectral trends 
to be seen. Sources with the strongest PAH features and weakest silicate absorption (sources 
15, 5, and 7) are among the lowest luminosity sources, in Figure 1. Sources of intermediate 
luminosity, in Figures 2 and 3, generally show conspicuous silicate absorption and weak 
PAH features. Finally, the highest luminosity sources (Figure 4) have the weakest silicate 
absorption and the weakest PAH features. 

As a comparison standard througho ut, we use the thoroughly studied ULIRG Markarian 
231 because its mid infrared spectrum (jWeedman et al.ll2005l : lArmus et al.l 120071 ) is similar 
in shape to many other absorbed sources. In Figure 5, we show this similarity by comparing 
the average spectrum of all the FSC sources in Table 2 having silicate absorption to the 
spectrum of Markarian 231. The average absorption depth is similar to that of Markarian 
231, for which the silicate optical depth of -0.7 corresponds to maximum extinction by the 
silicate feature of ~ 50% of the continuum flux. 

The distribution of continuum luminosities for the FSC sources is shown in Figure 
7, based on z/L,^(5.5 /im). This parameter is cho sen for compariso n to the large samples 
of ULIRGs, Syl, Sy2, and QSOs summarized in iHao et al.l (120071 ). All of these sources 
previously observed with the IRS have log[z/L;y(5.5 /xm)] < 46.1, and the median ULIRG 
luminosity is log[i^L,y(5.5 //m)] = 44.4. The luminosity limit of these previously observed 
so urces is exceeded by 3 of the FSC sources. The most luminous obs cured sources at z y 
2 jHouck et all I2OO5I : Ivan et al.ll2007h have log[z^L^(5.5 fim)] < 46.6 JPoUetta et al.ll2007h . 
and this extreme luminosity is matched by one of the FSC sources. The median luminosity 
of the FSC sources in Figure 7 is log[i/Li,(5.5 //m)] = 44.8. These results indicate, therefore, 
that the FSC sample encompasses a luminosity range which includes typical ULIRGs and 
extends to the most luminous sources discovered by Spitzer. 



3.3. Comparison of FSC Sample to Previous ULIRG samples 



One objective of the present study was to assemble a sample of IRAS sources substan- 
tially fainter in infrared and o ptical flux than the luminous ULIRG S deflned by th e IRAS 
1 Jy survey (IKim et al.lll998l) whose IRS spectra are discussed by llmanishi et al.l (120071 ). 
Farrah et all J2007h . and lPesai et al.l (120071 ). Objectives include learning how the nature of 



IRAS sources changes as flux becomes fainter, and especially to seek sources that are more 
similar in luminosity and redshift to the high-redshift sources discovered in Spitzer surveys 
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of sources at ~ 1 mJy. 



The FSC sample was chosen to have the faintest optical magnitudes that could be found 
in the FSC in order to select the IRAS sources with the most extreme values of observed 



infra red to optical fl ux ratio (IR/opt ). The faint Spitzer samples in Bootes (IHouck et al. 



20051 ) and the FLS (lYan et al.l 120071 ) have typical i? ~ 25 mag and /i/(25 /im) ~ 1 mJy. 
In the observer's frame, objects having similar IR/opt in the FSC, with typical /,y(25 fim) 
~ 100 mJy, should have i? > 20 mag. Very few FSC objects are so faint and so extreme 
in IR/opt, but all sources in Table 1 have photographic i? > 18. (We assume within the 
magnitude uncertainties that E and R are similar). 

Figure 8 compares the FSC sample and the ULIRG samples in 25 /im flux and redshift 
(using the Spitzer IRS fu{25 iim) for the FSC sources and the IRAS fu{25 //m) for the 
ULIRG sources). The FSC sample extends to higher redshifts than the previously published 
ULIRG samples which have IRS spectra. The maximum redshift of the FSC sample is 



0.93 w ith median redshif t of 0. 25 whereas the ULIRG samples tabulated by llmanishi et al. 



( 120071 ) and iFarrah et al.l (120071 ) have maximum z of 0.26. The combined dataset extends 
over a range of nearly 1000 in fiy{25 /im). 

The ratio of mid-infrared to optical or near-infrared flux is a measure of how much a 
source is dominated by dust absorption and dust continuum luminosity. In general, very 
dusty sources will have ratios enhanced both by increased extinction of shorter wavelengths 
and increased dust continuum at longer wavelengths, although detailed understanding of 
individual sources re quires modeling of the dust absorption and emission contributions (e.g. 
Marshall et al.l 120071 ). For comparing fluxes at different wavelengths among samples, homo- 
geneous all-sky data measuring optical magnitudes for IRAS sources are not available other 
than the photographic magnitudes used for our initial selection. However, the 2MASS survey 
( ISkrutskie et al.ll2006l ) encompasses all sources in the ULIRG samples and most sources in 
our FSC sample with near-infrared J,H,K photometry. 

In Figure 9, the flux distributions at 25 um and ratio /,y (25 pm ) / f,A J) for the FSC sampl e 
are compared with the ULIRG samples in llmanishi et al.l (120071 ) and iFarrah et al.l (120071 ). 
(The flux transformation adopted for 2MASS J mag is that zero magnitude corresponds to 
1590 Jy.) The median for the combined ULIRG samples is fu{25 fim)/ f^{J) = 2.42 and is 
the same for the FSC sample, counting limits. Although the FSC sample was chosen for faint 
optical magnitudes, this similarity in /i,(25 /im)//y(J) indicates that the FSC sample does 
not select in favor of sources with a greater dust content compared t o the ULIRG samples . 
The FSC sample does contain sources at greater redshifts than the llmanishi et al.l (120071 ) 
and IFarrah et al. (120071) ULIRG samp les, as shown in Fig. 8, although this result arises in 
part because the llmanishi et al.l (120071 ) sample was restricted to ULIRGs with z < 0.15. 
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3.4. Luminosities from AGN and Starbursts 



One major objective of previous ULIRG studies has been to determine the relative 
contributions of starbursts and AGN to the total luminosity of a source. A division of 
luminosity between starbursts and AGN at rest frame ~ 6 /um can be determined by the 
strength of the PAH features comp ared to the continuum, as measured by the equivalent 



width (EW) of the PAH features (IGenzel et al.l Il998l : iLaurent et all I2OOOI : llmanishi et al. 



20071 : iFarrah et al.l 120071 : iDesai et al.l 120071 ) . Pure st arbursts with no in dicators of AGN at 
any wavelength have rest frame EW(6.2 fim) > 0.45 (IBrandl et al.ll2006l ). Smaller equivalent 
widths correspond to an increasing continuum, assumed to be caused by an A GN. Very few 
U LIRGS have such s trong PAH features; only 2 of 49 in llmanishi et al.l (120071 ) and 8 of 107 



m 



Desai et al.l (120071 ). which indicates a substantial AGN contribution to the continuum at 



6 fiui . 

The distribution of PAH strength for the FSC sample is shown in Figure 7; only one of 
28 sources has 6.2 ^m EW > 0.45 and only three have EW > 0.2. These results indicate that, 
as in the ULIRG samples, the near-infrared continuum is generally much stronger relative to 
the PAH feature than in pure starbursts. It is also seen in Figure 7 that the most luminous 
sources at rest frame 5.5 /im in the FSC sample are sources with the smallest PAH equivalent 
widths. The simplest interpretation of these results is that the most luminous near-infrared 
sources are dominated by AGN power. 

How the luminosity divides in the near-infrared between AGN and starbursts is not 
necessarily the same as the separation of starburst and AGN power for bolometric luminosi- 
ties. Many efforts have been made to separate the bolometric contributions within ULIRGs 
of AGN and starbursts by assuming temp lates for the total AGN spectrum and the total 
starburst spectrum (e.g. iFarrah et al.ll2003l ). In general, results are that the far infrared lu- 
minosity from cooler dust is attributed primarily to the starburst whereas the near infrared 
luminosity from warmer dust is attributed to an AGN. The warmer AGN contribution leads 
to stronger near infrared continuum than from starbursts, causing the low equivalent widths 
for PAH features. 

Assumed templates can be incorrect if, for example, there is cooler dust outside of the 



cally thick clouds t hat produce hotter dust (e.g. 



Levenson et al. 


2007; 


Imanishi et al. 


2007 



PoUetta et al.l 120071 ). The templates may be examples of observed sources, but the intrinsic 
nature of the template sources may not be fully understood. It is desirable, therefore, to 
determine an estimate of AGN and starburst contributions which is independent of template 
assumptions. 
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For the rest frame parameters z/L,^(5.5 fim) and L(6.2 /im), which are measured for 
our sources in Table 3, empirical determinations relate these parameters to the bolomet- 
ric lum inosities of obscure d AGN [Ljr(AGN)] and of starbursts [Ljr(SB)]. For local star- 
bursts (IBrandl et al.l |2006| ) and high redshift submillimeter galaxies which are starbursts 
( Pope et al. 20081 ). an empirical relation has been found to be log[Lj,.(SB)] = log[L(6.2 jj,m)] + 
2.7±0.1; L(6.2 /xm) is the total luminosity of the 6.2 PAH feature, after subtracting the 
underlying continuum. 

This result derives from an observational comparison of PAH luminosity and bolometric 
luminosity; while there may be some extinction of the 6.2 fim feature, effects of extinction 
would be included in the empirical result. Adopting this relation also assumes that L(6.2 
yum) is measured the same in all sources. The uncertainty hsted does not include the ~ 10% 
systematic uncertainty in measurement of 6.2 um fluxes that m ay arise from differences in 
the fitting procedure discused i n section 2.2 . iPope et al.l (120081 ) state that their measures 
are "consistent" with those of iBrandl et al.l (120061 ). to which we also compare in section 
2. 2. Because our me asures agree within 10% of the Brandl measures, we adopt the relation 



m iPope et al.l (l2008l ) without applying a systematic correction. This relation applies to a 
"pure" starburst, with no AGN contribution to infrared luminosity. 

For AGN, it is assumed that the infrared continuum associated with the AGN arises 
from war m dust withi n a to rus surrounding the accretion regio n. Using the c l umpy torus 
model of iHonig et al.l (120061 ). empirical fits to overall SEDs in iPoUetta et al.l (120071 ) give 
log[Ljr(AGN)] = log[z^Li,(6.0 /um)] + 0.32±0.06. We make a nominal modification to this by 
transforming z/Li,(5.5 fim) to iyL,y{6.0 /jm) using Markarian 231 because of its similarity to 
the average spectrum of our absorbed sources in Figure 10. With this correction, we have 
log[Ljr(AGN)] = log[z/Li,(5.5 /um)] + 0.33±0.06. We adopt this relation as applying to a 
"pure" AGN, with no contribution to infrared luminosity from any source other than the 
torus heated by the AGN. 

These two relations for deriving Ljr(SB) and Ljr(AGN) are used to determine the star- 
burst and AGN luminosities for the FSC sources. Results are in Table 3 and Figure 10. We 
find that sources divide equally between those dominated by AGN luminosity and those by 
starburst luminosity. It is notable in Figure 10 that the median luminosity in the sample 
attributed to AGN, log Ljr(AGN) ~ 45.2, is the same as that attributed to starbursts, log 
Ljr(SB), although AGN sources reach to higher luminosities. 
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3.5. Radio Characteristics 



Because this F SC sample was selected with the requirement of a FIRST radio detection 
(IWhite et al.l 119971 ) . all sources have 1.4 GHz detections. We list in Table 2 the values of 
the parameter q, defined as q = log[/j,(25 fxm)/ f^{lA GHz)] in the observed frame. For 
determining q, sources are assumed to be unresolved so the peak /i/(1.4 GHz) is used. For 
the FSC sources, the peak flux typically differs by ~ 10% from integrated fluxes. As discussed 
in section 3.1, the much higher S/N and higher spatial resolution of the Spitzer detections 
gives more accurate photometry at 25 fim than the IRAS FSC fluxes so the IRS /;y(25 fim) 
are used to determine q. 

The comparison of q and z/Lj,(5.5 /xm) is shown in Figure 11. Larger values of q cor- 
respond to relatively weaker radio sources. The median value of q for the FSC sample is 
1.25. This is larger than the median o f q = 0.8 for faint Sp itzer First Look Survey sources 
detected at both 24yum and 1.4 GHz (lAppleton et al.ll2004l ). whi ch is the q value exp ected 



from previously known radio-infrared correlations for starbursts (ICondon et al.lll982l ). The 
median q is even less, q = 0.4, for faint sources in the Spitzer First Look Survey which have 
IRS spectra and z > 1 in IWeedman et al.l (j2006d ) . Some of these differences may be redshift 
effects. For example, Markarian 231 (shown in Figure 9) has q = 1.6 but, if at z = 2, would 
have q = 0.6. 

The results in Fig. 11 do not indicate that the value of q can be used as a discriminant 
between starburst and AGN sources. A median q of ~ 1.3 is found both for the starburst 
sources that have PAH features and for the remaining sources dominated by AGN. Also, 
the range of q is similar for both categories of sources. All of the values of q for the FSC 
sources are much la rger than in radio- loud AGN or quasars, for which q < -1 at any redshift 
(IHigdon et al.ll2005l ). The radio data, therefore, do not indicate any evidence of radio- loud 
AGN among the FSC sources. 



4. Summary 

The IRAS Faint Source Catalog (FSC) was used to select the 31 optically faintest FSC 
sources which are also identified in the FIRST radio survey, representing the optically faintest 
1% of IRAS extragalactic sources. 28 of these sources have been observed with the Infrared 
Spectrograph on Spitzer. Spitzer fluxes indicate that about half of the sources have fu{25 
fim) which are below the IRAS FSC flux limit. This indicates that about half of the optically 
faint sources taken from the FSC are actually spurious FSC detections. These sources are 
nevertheless detected by Spitzer because the presence of a FIRST radio source also provides 
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an infrared source of sufficient brightness for a Spitzer detection. 

This FSC sample reaches to higher redshifts and higher luminosities than the IRAS- 
discovered ULIRG samples observed spectroscopically with Spitzer, and overlaps the lumi- 
nosity range of 5'pii^er-discovered 24 //m sources at z ~ 2. The FSC sources have 0.12 < z 
< 1.0 and luminosities 43.3 < log[i/L,^(5.5 //m)] < 46.7. 15 of the sources have detectable 
PAH features, and 19 have measurable silicate absorption. Median properties of the sample 
having silicate absorption are very similar to the ULIRG Markarian 231 in silicate strength, 
continuum luminosity, ratio f^{25 iJ,m) / f^{J) , and relative radio luminosity. 

PAH luminosities are used to determine the starburst luminosity within each source, 
and predictions from dusty torus models are used to determine the AGN luminosity. Sources 
have similar bolometric luminosities arising from starbursts and from AGN and are equally 
divided between sources dominated by starbursts and sources dominated by AGN. The ratio 
of infrared to radio flux is not a measure of whether sources are dominated in the infrared 
by starburst or AGN luminosity. 

We thank Don Barry for technical assistance. This work is based on observations made 
with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, Cal- 
ifornia Institute of Technology under NASA contract 1407. Support for this work by the 
IRS GTO team at Cornell University was provided by NASA through Contract Number 
1257184 issued by JPL/Caltech. Support was also provided by the US Civihan Research 
and Development Foundation under grant ARP1-2849-YE-06. This research has made use 
of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under contract with the National Aeronautics 
and Space Administration. 
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Table 1. IRS Observations of Sources 



i>l UIIlUcI 


Source Name* 


CUUI QlIlclXoo 


A OR 


pi Ugl d-IIl 


c 

ximG 


rlQ<-eid 
Qaie 














s 


mo/day/yr 


1 


FSC 07247+6124 


072912.10+611853.5 


17540096 


30121 


120,60 


11/16/06 


2 


FSC 09105+4108 


091345.28+405632.6 


7771392 


1018 


120,180 


11/29/03 


3 


FSC 09121+2430 


091501.71+241812.2 


17540352 


30121 


84,60 


05/01/07 


4 


FSC 09235+5425 


092703.07+541206.6 


17540608 


30121 


120,60 


11/17/06 


5 


FSC 09284+0413 


093101.27+035955.2 


17540864 


30121 


84,60 


06/09/07 


6 


FSC 09425+1751 


094521.36+173753.4 


17541120 


30121 


84,60 


06/08/07 


7 


FSC 10219+2657 


102447.39+264209.0 


17541888 


30121 


84,60 


12/25/06 


8 


FSC 11257+5113 


112832.73+505721.1 


17542400 


30121 


120,60 


05/02/07 


9 


FSC 13080+3237 


131028.64+322049.2 


17542912 


30121 


120,60 


06/27/06 


10 


FSC 13297+4907 


133150.54+485150.6 


17543168 


30121 


120,60 


06/27/06 


11 


FSC 14448-0141 


144727.54-015330.3 


17543424 


30121 


120,60 


07/25/06 


12 


FSC 


14475+1418 


144954.86+140610.5 


17543680 


30121 


120,60 


07/29/06 


13 


FSC 


14481+4454 


144953.70+444150.3 


17543936 


30121 


120,60 


07/04/06 


14 


FSC 


14503+6006 


145135.04+595437.6 


17544192 


30121 


120,60 


07/04/06 


15 


FSC 


14516+3851 


145335.96+383913.1 


17544448 


30121 


120,60 


06/24/06 


16 


FSC 


14589+2329 


150113.19+232908.2 


4169216 


49 


56,120 


02/07/04 


17 


FSC 


15065+3852 


150825.42+384122.1 


17544704 


30121 


120,60 


07/01/06 


18 


FSC 


15307+3252 


153244.05+324246.7 


4983552 


105 


360,480 


03/04/04 


19 


FSC 


15385+4320 


154014.02+431042.4 


17544960 


30121 


120,60 


06/25/06 


20 


FSC 


15458+0041 


154823.38+003212.8 


17545216 


30121 


120,60 


03/20/07 


21 


FSC 


15492+3454 


155108.86+344533.6 


17545472 


30121 


120,60 


06/25/06 


22 


FSC 


15496+0331 


155206.16+032244.0 


17545728 


30121 


120,60 


03/16/07 


23 


FSC 


15585+4518 


160003.29+451046.2 


17545984 


30121 


120,60 


07/24/06 


24 


FSC 


16001+1652 


160222.38+164354.3 


17546240 


30121 


120,60 


03/10/07 


25 


FSC 


16073+0209 


160949.75+020130.8 


17546496 


30121 


120,60 


03/18/07 


26 


FSC 


16156+0146 


161809.36+013922.1 


17546752 


30121 


120,60 


09/17/06 


27 


FSC 


16242+2218 


162626.00+221145.9 


17547008 


30121 


120,60 


09/17/06 


28 


FSC 


17233+3712 


172507.40+370932.1 


4986880 


105 


360,480 


06/07/04 



'^Source name from the IRAS Faint Source Catalog. 

^Coordinates for the FIRST source which is closest to the FSC source. 

'^Total integration time for IRS short low modules (first number) and long low modules (second 
number) . 

'^Date of IRS observation 



Table 2. Observed Properties of Sources 



No. 


Source 




z 


f 1^(25 f^m) 


/„(25 Mm) 


E 






/„{1.4 GHz) 




/iy(5.5 Mm) 


/„(15 fim) 


EW{6.2fim) 


f(6.2fim)'^ 


f{11.3Mm)'' 












IRAS mjy 


IRS mJy 


mag 


mag 


mJy 




m Jy 


m Jy 










1 


07247 


0.137 


175 


119 


18. 





16.08 


3.13 


1.57 


4.1 


47.2 


0.08±0.01 


2.94±0.42 


2.48±0.46 


-0.31 


2 


09105 


0.446 


333 


333 


18. 


1 


16.83 


5.79 


1.76 


54 


297 


< 0.002 


< 1.3 


< 0.2 


-0.34 


3 


09121 


0.84^ 


179 


100 


20. 





16.56 


9.80 


1.01 


42.5 


119 


< 0.01 


< 1.8 


< 0.4 


0.0 


4 


09235 


0, 


,123 


136 


97 


18. 


3 


16.75 


3.31 


1.47 


2.0 


42.8 


0.19±0.02 


4.62±0.47 


2.67±0.16 


0.0 


5 


09284 


0, 


,146 


147 


25 


18. 


2 


16.11 


1.05 


1.38 


1.5 


6.9 


0.56±0.03 


12.14±0.73 


10.0±0.2 


0.0 


6 


09425 


0, 


,130 


436 


322 


18. 


6 


15.72 


38.7 


0.92 


17.8 


162 


< 0.1 


< 0.3 


< 3.91 


-0.78 


7 


10219 


0, 


,225 


151 


28 


18. 


8 


16.77 


1.68 


1.22 


1.3 


11.8 


0.29±0.04 


4.53±0.7 


3.92±0.5 


0.0 


8 


11257 


0, 


,197 


123 


75 


18. 


4 


16.55 


1.77 


1.63 


3.0 


35.8 


0.07±0.01 


3.12±0.6 


1.91±0.4 


-0.51 


9 


13080 


0, 


,997"" 


271 


9.0 


19. 


2 


16.34 






3.5 


10.3 








0.0 


10 


13297 


0, 


,128 


139 


102 


18. 


5 


16.68 


1.17 


1.94 


3.6 


53.1 


< 0.05 


< 1.8 


< 1.5 


-0.69 


11 


14448 


0, 


,210 


176 


38 


18. 


,9 


> 17.7 


1.78 


1.33 


3.8 


23.4 


< 0.02 


<0.9 


2.01±0.65 


-0.62 


12 


14475 


0, 


,251 


143 


114 


18. 


5 


16.68 


12.8 


0.95 


11.6 


91 


0.03±0.005 


3.4±0.52 


< 3.2 


-0.89 


13 


14481 


0, 


,670 


85 


64 


21. 


,0 


> 17.7 


10.7 


0.78 


19.2 


67.6 


< 0.006 


< 1.1 


< 0.5 


0.0 


14 


14503 


0, 


,577 


79 


58 


21. 


,0 


> 17.7 


16.3 


0.55 


8.4 


60 


< 0.02 


< 1.9 


< 1.0 


-1.01 


15 


14516 


0, 


,153 


85 


24 


19. 


,3 


16.77 


2.07 


1.06 


0.75 


9.4 


0.20±0.03 


2.24±0.38 


2.49±0.14 


-0.87 


16 


14589 


0, 


,261 


84 


53 


17. 


,8 


15.98 


4.49 


1.07 


11.7 


37.8 


< 0.017 


< 1.73 


< 3.4 


-0.24 


17 


15065 


0, 


,355 


82 


45 


20. 


2 


> 17.7 


2.43 


1.27 


4.9 


32.9 


< 0.013 


< 0.84 


4.41±0.27 


-0.59 


18 


15307 


0, 


,927 


71 


46 


19. 


,0 


> 17.7 


5.7 


0.91 


9.4 


70.0 


< 0.006 


< 0.5 


< 2.4 


-0.32 


19 


15385 


0, 


,380 


92 


76 


18. 


5 


> 17.7 


1.10 


1.84 


6.0 


54.6 


0.02±0.005 


1.22±0.34 


< 1.1 


0.0 


20 


15458 


0, 


,254 


126 


40 


18, 


7 


16.70 


2.53 


1.20 


4.3 


25.2 


< 0.01 


< 0.59 


< 0.26 


0.0 


21 


15492 


0, 


,311 


80 


40 


19. 


,6 


> 17.7 


1.08 


1.57 


6.0 


29.7 


< 0.02 


< 0.9 


1.56±0.42 


-1.46 


22 


15496 


0, 


,193 


144 


122 


18. 


5 


16.20 


28.5 


0.63 


13.7 


86.7 


< 0.01 


< 1.8 


2.39±0.47 


-1.10 


23 


15585 


0, 


,486 


64 


11.0 


20. 


,0 


> 17.7 


1.58 


0.87 


3.1 


11.9 


0.10±0.04 


2.0±0.82 


3.25±0.51 


-2.2 


24 


16001 


0, 


,672 


144 


104 


18. 


1 


16.62 


1.90 


1.74 


22.9 


109 


< 0.01 


< 0.25 


< 1.29 


0.0 


25 


16073 


0, 


,223 


153 


123 


19. 


2 


17.46 


5.41 


1.36 


9.6 


96 


0.06±0.002 


7.5±0.3 


< 4.1 


-1.5 


26 


16156 


0, 


,133 


279 


261 


18. 


,3 


16.57 


7.9 


1.52 


27 


154 


< 0.01 


< 2.8 


< 0.8 


-2.6 


27 


16242 


0, 


,157 


71 


31 


19. 


5 


> 17.7 


1.39 


1.35 


1.3 


15.1 


< 0.01 


< 0.25 


1.07±0.12 


-2.7 


28 


17233 


0, 


702 


50 


29 


21, 


,0 


> 17.7 


3.16 


0.96 


1.7 


34.3 


< 0.04 


< 0.78 


< 1.2 


-0.51 



^Photographic E magnitude from APS. 

magnitude from 2MASS (Skrutskic ct al. 2006). 
'^q — log[/jy(25 /J,m)//jy (1.4 GHz)] in observed frame, using /;y(l.4 GHz) of "peak flux" from FIRST, corresponding to an unresolved source. 

'^Total flux of feature in units of 10 ^"'"W cm ^, fit with single Gaussian. Uncertainties in fluxes and equivalent widths are statistical uncertainties from fitting of feature and do 
not include systematic uncertainties of < 10% which may arise from fitting procedure chosen, as discussed in text. 

^Tgi is measure of depth of 9.7 ^m absorption feature, defined as T^i — ln[/j^ (abs)//jy (cont)] , for /j^ (abs) measured at the maximum depth of absorption and /jy(cont) the 
unabsorbed continuum as extrapolated in lSpoon et aT ] {2003). 

^Redshift for no. 3, a red quasar, is optical redshift from fUrrutia et al.| j2005|) . and for number 9, a blazar, is from fPerry et al.| j 19781 ) ; other redshifts are determined from the 
IRS spectra. 



Table 3. Luminosities of Sources 



Number 


Source 


Mpc 


ergs s ^ 


log[[vLiy(15 ^m}] 
ergs s ^ 


log[L(fa.2)J 
ergs s 


ergs s ^ 


log[L(bB)J 
ergs s 


^og[ij(,AtjiN )\ 
ergs s ^ 




1 


07247 


637 


43.98 


44.60 


42.15 


42.42 


44.85 


44.31 


> 0.78 


2 


09105 


2460 


46.17 


46.47 


<42.96 


<43.03 


<45.66 


46.50 


< 0.13 


3 


09121 


5340 


46.63 


46.64 


<43.79 


<43.88 


<46.49 


46.96 


< 0.25 


4 


09235 


565 


43.58 


44.46 


42.25 


42.44 


44.95 


43.91 


> 0.92 


5 


09284 


683 


43.60 


43.83 


42.83 


43.09 


45.53 


43.93 


> 0.98 


6 


09425 


602 


44.57 


45.09 


<41.10 


<42.26 


<43.80 


44.90 


< 0.07 


7 


10219 


1108 


43.93 


44.45 


42.82 


43.09 


45.52 


44.26 


> 0.95 


s 


11257 


952 


44.17 


44.81 


42.53 


42.73 


45.23 


44.50 


> 0.84 


9 


13080 


6610 


45.70 


45.73 












10 


13297 


594 


43.87 


44.60 


<41.87 


<42.14 


<44.57 


44.20 


< 0.70 


11 


14448 


1025 


44.34 


44.68 


<42.05 


<42.56 


<44.75 


44.67 


< 0.55 


12 


14475 


1254 


44.98 


45.43 


42.80 


<43.09 


45.50 


45.31 


> 0.61 


13 


14481 


4036 


46.08 


46.20 


<43.33 


<43.50 


<46.03 


46.41 


< 0.29 


14 


14503 


3360 


45.59 


46.01 


<43.41 


<43.59 


<46.11 


45.92 


< 0.60 


15 


14516 


719 


43.34 


44.00 


42.14 


42.46 


44.84 


43.67 


> 0.94 


16 


14589 


1310 


45.01 


45.09 


<42.55 


<43.02 


<45.25 


45.34 


< 0.44 


17 


15065 


1870 


44.92 


45.31 


<42.54 


43.34 


45.24 


<45.25 


< 0.50 


18 


15307 


6030 


46.06 


46.50 


<43.32 


<44.09 


<46.02 


46.39 


< 0.30 


19 


15385 


2030 


45.07 


45.59 


42.78 


<43.06 


45.48 


45.40 


> 0.55 


20 


15458 


1270 


44.55 


44.89 


<42.05 


<42.21 


<44.75 


44.88 


< 0.43 


21 


15492 


1606 


44.89 


45.14 


<42.43 


<42.87 


<45.13 


45.22 


< 0.45 


22 


15496 


930 


44.81 


45.18 


<42.26 


<42.63 


<44.96 


45.14 


< 0.40 


23 


15585 


2730 


45.01 


45.15 


43.25 


43.67 


45.95 


45.34 


> 0.80 


24 


16001 


4050 


46.17 


46.41 


<42.69 


<43.48 


<45.39 


46.50 


< 0.07 


25 


16073 


1090 


44.78 


45.35 


43.03 


<43.22 


45.73 


45.11 


> 0.80 


26 


16156 


620 


44.78 


45.09 


<42.11 


<42.22 


<44.81 


45.11 


< 0.34 


27 


16242 


740 


43.59 


44.24 


<41.22 


<41.94 


<43.92 


43.92 


< 0.50 


28 


17233 


4270 


45.07 


45.94 


<43.23 


<43.64 


<45.93 


45.40 


< 0.77 



^^uminosity distance determined by E.L. Wright, ^^http: //www. astro. ucla.edu'''^wright/Cosn:ioCalc.htmT| for Hq — 71 km s "'" Mpc "'" , f2j\/— 0.27 and f2^— 0.73. 
•^Total luminosity of starburst, assuming that log[Li^(SB)] = log[L(6.2)]+ 2.7±0.1, from Brandl et al.{2006) and Pope et al.(2007). 

■^Total luminosity from dusty torus of AGN, assuming log[L^^ (AGN)] — log[fLi^ (5.5 /^m)] + 0.33ib0.06, from Pollctta et al.(2007) and assuming that all continuum 
luminosity i/'Li^(5.5 ixra) arises from an AGN. 

■^Fraction of total luminosity which arises from a starburst. If a PAH feature is detected, this fraction is shown as a lower limit because some of the continuum 
at 5.5 /Am which is attributed to an AGN actually arises from the starburst. If no PAH feature is detected, this fraction is shown as an upper limit because there is 
no direct evidence that a starburst is present. Source 9 is a blazar so the parameters for determining Lj,,(AGN) and L^^(SB) are not applicable. 
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Fig. 1. — IRS spectra of FSC sources in rest frame, ordered by i/L^{5.5 /im) in Table 3, 
lowest luminosity at the top. Sources of higher luminosity continue in Figure 2. 
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Fig. 2. — IRS spectra of FSC sources in rest frame, ordered by z/L^(5.5 /im) in Table 3, 
lowest luminosity at the top, continued from Figure 1. Sources of higher luminosity continue 
in Figure 3. 
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Fig. 3. — IRS spectra of FSC sources in rest frame, ordered by i/L^{5.5 /xm) in Table 3, 
lowest luminosity at the top, continued from Figure 2. Sources of higher luminosity continue 
in Figure 4. 
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Fig. 4. — IRS spectra of FSC sources in rest frame, ordered by z/L,^(5.5 /xm) in Table 3, 
lowest luminosity at the top, continued from Figure 3. Source having highest luminosity is 
number 3. 
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Fig. 5. — Comparison of Markarian 231 (top) with average of the 19 FSC galaxies in Table 
1 with measured silicate absorption (bottom), arbitrarily normalized. 
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Fig. 6. — Comparison of /,y(25 /iin) measured with Spitzer IRS to /jy(25 yum) listed in IRAS 
FSC for sources in Table 2. Uncertainties in IRS fluxes arise primarily from calibration 
uncertainties and are ~ ± 5%, comparable to the size of plotted symbols. Large discrepancies 
between IRS and IRAS fluxes for fainter sources are discussed in text as indication that many 
IRAS sources are spurious. 
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Fig. 7. — Comparison of luminosities and PAH strengths for FSC sources; filled squares are 
sources with PAH detections and open squares show PAH upper limits for sources without 
PAH detections. The star indicates Markarian 231. Error bars show uncertainties in PAH 
EWs from Table 2. 
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Fig. 8. — Distribution of redshift and /jy(25 fim) for ULIRG samples in Farrah et al. 
(2007) (filled triangles) and Imanislii et al. (2007) (filled circles) compared to FSC sam- 
ple in present paper(filled squares). The star indicates Markarian 231. The fu{25 fim) 
shown for ULIRGS are taken from IRAS fluxes; the fu{25 fim) shown for FSC sources are 
synthetic /,y(25 /im) that should be observed by IRAS as determined from Spitzer IRS spec- 
tra. Uncertainties in IRS fluxes arise primarily from calibration uncertainties and are ~ ± 
5%, comparable to the size of plotted symbols. 
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Fig. 9. — Comparison of mid- infrared and near-infrared (2MASS J) fluxes for ULIRG samples 
in Imanislii et al. (filled circles), Farrali et al. (filled triangles), and FSC sample in Table 2 
(squares); open squares are limits for sources without J detections so the f(25 ^m)/f(J) ratio 
is greater than shown. The star indicates Markarian 231. The fiy{25 iim) shown for ULIRGS 
are taken from IRAS fluxes; the fu{25 /im) shown for FSC sources are synthetic /i/(25 ^m) 
that should be observed by IRAS as determined from Spitzer IRS spectra. Uncertainties in 
IRS fluxes arise primarily from calibration uncertainties and are ~ ± 5%, comparable to the 
size of plotted symbols. 
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Fig. 10. — Comparison of bolometric luminosities from AGN with those from starbursts for 
FSC sample in Table 3, assuming that log[Ljj.(SB)] = log[L(6.2 /im)]+ 2.7 for PAH lumi- 
nosity L(6.2) and that log[Ljr(AGN)] = log[z^L,^(5.5 //m)] + 0.33 for continuum luminosity 
i^Lj,(5.5 /im). Filled squares are starbursts with PAH detections so the L(AGN) derived from 
continuum luminosities are upper limits for L(AGN) because some continuum may arise from 
the starburst; open squares are sources without PAH detections so the L(SB) derived from 
PAH luminosities are upper limits for L(SB) determined from upper limits on L(6.2). The 
star indicates Markarian 231. 
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Fig. 11. — Comparison of luminosities and radio strength for FSC sources; q = log[/,y(25 
/im)//j,(1.4 GHz)] in observed frame. Filled squares are starburst sources (as defined by PAH 
detections), and open squares are AGN (no PAH detections). The star indicates Markarian 
231. Uncertainties in q are typically ± 10%, shown by the representative error bar, arising 
primarily from uncertainty whether to adopt "peak" flux or "integrated" flux for the FIRST 
sources. 



